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S.J. Lee 1 ' 2 and A.Z. Mekjian 2 
1 Department of Applied Physics, Kyung Hee University, Yongin-si, Gyeonggi-do 446-701, Korea and 
2 Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08854, USA 

The nuclear binding energy is studied using a finite temperature density functional theory. A 
Skyrme interaction is used in this work. Volume, surface, and symmetry energy contributions to 
the binding energy are investigated. The case of neutron skin is considered in detail. The neutron 
skin modifies the mass A dependence of various terms and / dependence of the skin thickness is 
proportional to / for the case of same central density. 



I. INTRODUCTION 



The mass formulae [H-Q characterizes the binding energy of nuclei in terms of the proton number Z and neutron 
number N or nucleon number A = Z + N. Volume, surface, Coulomb and pairing terms appear. Of importance and 
one of the least determined parts of the binding energy expression is the symmetry energy term. The symmetry energy 
has both a volume symmetry part and a surface symmetry part, similar to the division of the volume and surface 
terms. The division into volume and surface energies occurs because observed nuclei have finite nucleon number with 
■ a maximum value of A ~ 250. The symmetry energy appears in the total binding energy with factors involving both 
(f) ' ((N — Z)/A) 2 A and ((N — Z)/A) 2 A 2 / 3 for the volume and surface symmetry terms respectively. For a heavy nucleus 
with A = 200, Z = 80, ./V = 120 the (N - Z)/A = 1/5 and {(N - Z)/A) 2 = 1/25. The isospin dependence of the 
binding energy arises in part from kinetic energy differences between protons and neutrons and from interaction terms 
arising from the isospin dependence of the nuclear force between nucleons. This interaction has terms involving the 
fi ■ Tj which arises from the exchange of isovector mesons. The surface terms in the mass formulae arise from the loss 
of binding energy for nucleons near the surface of a nucleus. In a neutron star, the volume term only exists. Thus, 
of importance for neutron star physics is the extraction of the volume term from properties of known finite nuclei 
which contain both volume and surface terms. The symmetry term also plays a significant role in heavy ion collisions 
043 1 m neutron star physics in the valley of nuclear stability and associated neutron and proton drip lines, in 
neutron halo nuclei in isobaric analog states [9l-fllj|. and in giant isovector dipole states Q. Extensive discussion of 
the symmetry energy can be found in the work of Danielewicz et al. fl2| . The role of the symmetry energy in the 
nuclear surface has also been studied in Ref . [la . An extensive review can be found in Ref. M . 
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p | In our preveous paper [13| , we have examined the T dependence of the expansion coefficients in mass formular by 

I/-) ■ minimizing the Hclmholtz free energy using Skyrme interaction. In there we have assumed R n = R p and thus missed 
any effects of neutron skin which would exists in asymmetric nuclei. We estimated approximately neutron skin effect 
for T 2 term in kinetic energy by setting the central densities to be the same p nc — p pc and found the effect of R n ^ R p 
' is small. But this is only a part of neutron skin effect thus we will study the neutron skin dependence more fully here 
t-H . by expanding energy with the neutron skin size R n — R p . 

In this paper we study the neutron skin effect on the volume and surface contributions to the symmetry energy. Our 
. y—{ \ approach is based on a finite temperature density functional theory and we use a Skyrme type of interaction which 
^ ■ we develop in the next section. The same Skyrme approach was also used to study phase transition in Ref. (l4l - fl7t . 
The symmetry energy and the neutron skin effects on symmetry energy and other energy coefficients for three cases 
of Skyrme interactions are studied in Sect. IIIII and the results are summarized in Sect. IIV1 



II. BINDING ENERGY IN A DENSITY FUNCTIONAL APPROACH. 



A density functional theory based on a Skyrme interaction will be used in our investigation which is summarized 
in Appendix [Al The density distribution is taken to be of the Saxon- Wood form: 

Pq (r) = p q (R q ) = — -ggL_ (1) 

The q = p,n for protons, neutrons. We will allow the central density p qc and the R q for protons and neutrons to be 
different in general. The diffuseness a q are taken to be the same. Two limiting cases can also be considered. These 
are: 1. R p = R n so that the difference in N ^ Z nuclei is in the central density or p pc ^ p nc and 2. R p ^ R n and 
Ppc — Pnc- In case 2 the neutron distribution reflects a neutron halo for N > Z . 

Using this density functional approach the binding energy of nuclei as a function of mass number A, proton number 
Z and temperature T can be evaluated. The Weizacker semiemperical mass formulae [l[, its extension by Myers and 
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Swiatecki [18-21] and also including finite temperature effects [13[ is 

E{A,Z,T) = E V {T)A + E S {T)A 2 ' 3 + S V {T)I 2 A + S S {T)I 2 A 2 ' 3 

z 2 z 2 z 4 / 3 
+ E cjm + Edi f~A + e °*aP* + cAA ~ 1/2 (2) 

where / = (JV — Z)/A = (A — 2Z)/A. The B = —Ey is the usual bulk energy per nucleon. The Edif and E ex are 
the coefficients for the diffuseness correction and the exchange correction to the Coulomb energy. For the pairing 
correction with constant A, c = +1 for odd-odd nuclei, for odd-even nuclei, and —1 for even-even nuclei. The above 
formula at T = is the well known Weizacker semiempirical mass formula [H-Q studied extensively by Myers and 
Swiatecki p^ - [2~l| . Early studies excluded the surface symmetry term Ss and only the surface term E$ was included. 
The values of the coefficients as found in textbooks such as Ref. @,H are E v {0) = -B(0) ss -16, E s (0) « 17, and 
Sy(0) s» 24 in MeV. The ratio E$/B of surface to bulk energy at T — is very close to unity. 

In our preveous paper (l3j . we have examined the T dependence of the expansion coefficients in Eq.Q by minimizing 
the Hclmholtz free energy using a Skyrme interaction with the density distribution of Eq.([T]). With R n = R p = R, we 
can integrate Skyrme interaction analytically to obtain total energy as a function of R. Using this energy function, 
the total energy E(A, Z, T) minimizing free energy is found by varying the value of R for a nucleus with Z protons and 
iV neutrons at temprerature T. Then we use Eq.([2]) for various nuclei to obtain the expansion coefficients. Since we 
have assumed R n = R p in our previous studies, we missed any effects of neutron skin which would exist in asymmetric 
nuclei. We estimated approximately neutron skin effect for T 2 term in kinetic energy, Eq. (|A3[) . by setting the central 
densities to be the same, p nc = p pc — p c /2 where p c is the total central density obtained with R n = R p , and found 
the effect of R n ^ R p in kinetic energy is small. 

In the present paper we now consider neutron skin effects more fully. Specifically, we study the neutron skin size 
t = R n — R p dependence of total energy E(A, Z, T, t). Then we expand each coefficient Ei of Eq.J2|) as 

Ei(T,t) = Ei(T)+Ei(T) Bk -^ (3) 

Since we can integrate a Skyrme interaction analytically for the case of R n = R p = R, we can obtain the neutron skin 
t dependence by expanding the integral for the case of R n ^ R p around R n — R p — R. This can be done by expanding 
the density distribution of Eq.([T]) around R which we now discuss. It should also be noted that in Eq.([3]) the surface 
correction appears with the factor t/A 1 / 3 which originates from the dimensionless factor t/R with R = tqA 1 / 3 . The 
extra I /A 1 / 3 introduces a different A dependence from the non skin part Ei(T). 

From Eqs.Q and (J3j> we see that the surface thickness factor t/A 1 / 3 changes the A dependences in the mass 
formular. In particular AE l (T) sk t/ A 1 / 3 acts as a surface term as A 2 / 3 and A 2 / 3 E i {T) sk t/A 1 / 3 goes like A 1 / 3 when 
we look at the whole A dependences of the neutron skin size dependent coefficients. When R q w R 7 the neutron skin 
size t/R is proportional to -£w as can be seen in Appendix [B] Thus the skin dependence of the energy expansion 
coefficients have this extra factor of A dependence. On the other hand, when neutron and proton central densities 
are the same, the neutron skin size t/R is proportional to (N — Z)/A (see Appendix |B|) . Then the skin dependence 
of the energy expansion coefficients introduce an extra I = (N — Z)/A behavior and the energy expansion of Eq.Q 
with Eq.J3]) becomes a third power expansion in the isospin factor /. That is t/A 1 / 3 ~ I since t/R ^ I. For such a 
case, by expanding the empirical nuclear energy for various nuclei together with I and I 3 terms included, we may be 
able to extract the information about the neutron skin size from the odd term in I. 

The central density p qc of the density distribution Eq.([T]) for a given value of R q should be determined to give a 
fixed number of nucleons N q . This normaization condition gives the expansion of p qc (R q ) as given by Eq. (|B4|) in 
Appendix IB1 up to first order in t q = (R q — R). Thus the first order correction coefficient to /?™ due to fixed N q is 

Furthermore if we keep the particle number N q and the total central density p c = p nc + p pc to be a constant while 
varying R qi that is, 



the zcroth order term times 



This result is independent of which type of particle, neutron or proton. 



p c (R) = p nc {R„) + Ppc(Rp) = Pnc{R) + Ppc(R) (4) 

then we have 

Rn R ~P^' 

N 

t p = Rp — R = — —t (5) 
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to lowest order in t — R n — R p when we expand about R q = R (see Appendix |Bj . 

The Fermi density p q {r) is then expanded about R q = R up to first order in {R q — R) as given by Eqs. (|B18l) - (|B20[) . 
The results of Appendix [B] show that the quantity p™ (r) + p™ (r) has a first order correction from skin size t and the 
first order correction vanishes for m — 1. That is the total density p(r) — p n (r) + p p (r) is independent of the neutron 
skin size t = R n — R p up to first order. Thus only the explicit p q dependent terms in Skyrme interaction, not the 
total density p dependent terms, depend on the skin size t up to the first order. 

Since F[p q ] = J d 3 rf(p q (r)) for a Fermi density p q (r) = 1+e (, P -a g )/ n ; where f{p q ) is only a function of single density 
p q , then can be integrate exactly. The result is a function F(R q ) of R q [l3| . which we can expand easily in 

terms of t q = R q — R around F(R). That is 

F(R q ) - F(R) + -^- tq + + ■ ■ (6) 

This procedure is much simpler than using a method of expanding the density in t first and then integrating the 
results. For fixed N q , the central density p qc (R q ) is also a function of R q as in Eq. (|B4l) . The expansion of integral of 
various power of p q are summarized in Appendix [Cj 

For the integral of the form of J d 3 rp a p q n , we need to expand p q n {R q ) around R first then integrate each term which 
is now a function of R only. The results for various cases are also given in Appendix [Cj 



III. NETRON SKIN SIZE DEPENDENCE 



Here we examine the neutron skin size dependence of nuclear energy using various Skyrme interactions. We used 
three sets of Skyrme parameters with different values of the effective mass in symmetric nuclear matter which are 
SLy4 with m*/m = 0.69, SkM* with m*/m = 0.79, and SkM(m* = m) with m*/m = 1. The results for these three 
cases are given in Tableland summarized in the following equations. The three cases cover a wide range of interaction 
types in terms of effective mass within the many various Skyrme interactions. 

For Coulomb energy, from Eqs.dUg]) and dUS]) with a = 0.53 fm and R = 1.25A 1 / 3 fm, E c part is 0.6912000, E ex 
part is -0.5278064 and Ediff -1.430810 and the skin dependence of Coulomb energy Ec part is -0.5529600, E ex part 
is 0.4222251 and E dlff 3.433944, that is 



, , z 2 , z 2 , z 4 / 3 

E C (T, t) + E dlf (T, t) — + E ex (T, t) 



z 2 z 2 z^i 3 

0.6912— - - 1.4308— - 0.5278- 



A 1 / 3 ■ A ' A 1 / 3 



Z 2 Z 2 z 4 / 3 

-0.5530--7Tr + 3.4339— + 0.4222- 



A 1 / 3 A A 1 / 3 



(7) 



A 1 / 3 



in MeV and fm units. 

In TableU the items labeled "Ei(Tf and ll E l {T) sk " are the energy expansion coefficients given in Eqs.Q and ([3]). 
The items labeled by "T-indp", "T 2 ", and "Kine" under Ei(T) are the values of temperature T independent part, 
T 2 dependent term in kinetic energy (Eq. (|A3p ). and the total kinetic energy contribution to the energy expansion 
coefficients Ei(T) respectively. Similarly the items labeled by "T-indp", "T 2 " , and "Kine" under Ei(T) s k are the 
values of temperature T independent part, T 2 dependent term in kinetic energy fEq. (IA3[ )). and the total kinetic 
energy contribution to the neutron skin dependent part of the energy expansion coefficients Ei(T) s k respectively. 

The results of Table U show the following features for the skin size dependence of various coefficients Ei(T) s k 
of Eq.([3]) in the mass formulae which behave as Ei(T) s i c t/A 1 / 3 . One overall feature for all components Ei(T) s k is 
the weak dependence on temperature. The volume binding energy coefficient Ei{T) s k — Ev{T) s k of about 5^7 
MeV/fm has the smallest value of all the Ei(T) s k terms. Comparing to this value, the skin independent volume 
energy coefficient E\r{T) is about —15 MeV. The surface energy coefficient Es(T) s k has a somewhat larger values of 
about — 17 ~ —22 MeV/fm and is negative compared to about 20 MeV for skin independent coefficient Es(T). The 
neutron skin size dependence of the symmetry energy terms have the following features. The volume symmetry energy 
coefficient S\/(T) s k ~ 60 MeV/fm while the surface symmetry energy coefficient has the largest magnitude of about 
—600 ~ —700 MeV/fm. Comparing to these the neutron skin size independent coefficients are «SV(T) « 20 ~ 30 MeV 
and Ss(T) w 30 ~ 50 MeV. Some dependence on the effective mass is present for all coefficients as can be seen in 
comparing SLy4 (m*/m = 0.7), SkM* (m*/m = 0.8), and SkM(m* = m). However the dependences of the individual 
coefficients on the temperatures and on the Skyrme parameters we used are not so sensitive. The magnitude of 
the neutron skin size dependent coefficients are largely different between different coefficients ranging from about 
5 MeV/fm for volume energy coefficient to about 700 MeV/fm for surface symmetry energy coefficient while the 
magnitudes of the skin size independent parts were of the same order of magnitude ranging from about 15 MeV for 
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TABLE I: Neutron Skin Dependence of Energy coefficient minimizing free energy. 



SLy4 



SkM* 



SkM(m* = m) 















-15.308 
-15.308 
0.01151 
28.653 
5.5202 
5.5202 
-0.00722 
-3.8780 



-15.296 
-15.308 
0.01145 
28.594 
5.5100 
5.5172 
-0.00722 
-3.8654 



-15.263 
-15.308 
0.01121 
28.423 
5.4792 
5.5081 
-0.00723 
-3.8269 



-15.217 
-15.312 
0.01061 

28.124 
5.42745 

5.4927 
-0.00725 
-3.7576 



-15.127 
-15.127 
0.01925 
25.483 
5.5837 
5.5837 
-0.00761 
-0.4893 



-15.108 
-15.127 
0.01921 
25.444 
5.5730 
5.5806 
-0.00762 
-0.4995 



-15.051 
-15.127 
0.01908 
25.319 
5.5411 
5.5717 
-0.00763 
-0.5302 



-14.962 
-15.130 
0.01868 

25.105 
5.48746 

5.5563 
-0.00765 
-0.5816 



-15.310 
-15.310 
0.03948 
19.957 
6.9809 
6.9809 
-0.00793 
-0.3867 



-15.270 
-15.310 
0.03948 
19.978 
6.9701 
6.9780 
-0.00793 
-0.3963 



-15.152 
-15.310 
0.03950 
20.037 
6.9379 
6.9697 
-0.00793 
-0.4253 



-14.954 
-15.310 
0.03948 

20.134 
6.88390 

6.9553 
-0.00794 
-0.4737 



20.008 
20.008 
0.5483 
-32.106 
-17.077 
-17.077 
0.01959 
12.325 



20.559 
20.009 
0.5501 
-31.638 
-17.040 
-17.059 
0.01958 
12.290 



22.247 
20.024 
0.5559 
-30.255 
-16.928 
-17.006 
0.01957 
12.179 



25.204 
20.099 
0.5672 
-27.943 
-16.739 
-16.915 
0.01955 
11.981 



18.756 
18.756 
0.5371 
-26.275 
-17.261 
-17.261 
0.02082 
1.4093 



19.296 
18.757 
0.5388 
-25.810 
-17.222 
-17.243 
0.02082 
1.4377 



20.948 
18.771 
0.5443 
-24.402 
17.103 
-17.187 
0.02082 
1.5238 



23.836 
18.841 
0.5550 
-22.060 
-16.903 
-17.090 
0.02080 
1.6681 



18.303 
18.303 
0.4995 
-15.845 
-21.691 
-21.691 
0.02201 
1.1228 



18.804 
18.304 
0.5003 
-15.379 
-21.656 
-21.678 
0.02201 
1.1493 



20.321 
18.310 
0.5027 
-13.966 
-21.550 
-21.638 
0.02200 
1.2296 



22.904 
18.339 
0.5072 
-11.598 
-21.372 
-21.569 
0.02199 
1.3638 



31.113 
31.113 
0.5445 

-24.674 
63.203 
63.203 

0.03095 
1.3320 



31.664 
31.115 
0.5489 

-24.568 
63.225 
63.194 

0.03104 
1.3578 



33.401 
31.150 
0.5626 

-24.103 
63.266 
63.140 

0.03134 
1.3656 



36.645 
31.333 
0.5902 

-23.455 
63.306 
63.019 

0.03188 
1.3511 



29.655 
29.655 
0.5482 
-33.636 
62.787 
62.787 
0.02954 
-1.3400 



30.210 
29.657 
0.5527 
-33.481 
62.836 
62.806 
0.02962 
-1.2880 



31.960 
29.693 
0.5667 
-32.837 
62.956 
62.836 
0.02988 
-1.1223 



35.221 
29.875 
0.5940 
-31.703 
63.125 
62.852 
0.03037 
-0.8333 



19.685 
19.685 
0.4227 
-6.8813 
63.366 
63.366 
0.03144 
-0.8709 



20.110 
19.686 
0.4241 
-6.5197 
63.381 
63.350 
0.03151 
-0.8224 



21.412 
19.698 
0.4287 
-5.4723 
63.428 
63.301 
0.03171 
-0.6789 



23.695 
19.756 
0.4377 
-3.7475 
63.498 
63.209 
0.03207 
-0.4350 



-41.035 
-41.035 
-2.5769 
171.020 
-711.559 
-711.559 
-0.91606 
-276.281 



-43.642 
-41.046 
-2.5962 
170.079 
-711.507 
-710.590 
-0.91767 
-275.271 



-51.815 
-41.196 
-2.6548 
165.833 
-711.128 
-707.436 
-0.92289 
-271.698 



-66.958 
-41.984 
-2.7749 
159.697 
-710.053 
-701.662 
-0.93232 
-265.302 



-43.596 
-43.596 
-2.5598 
166.818 
-706.246 
-706.246 
-0.94030 
-40.3526 



-46.185 
-43.606 
-2.5793 
165.864 
-706.252 
-705.310 
-0.94192 
-41.5910 



-54.328 
-43.767 
-2.6404 
162.452 
-706.059 
-702.271 
-0.94708 
-45.3514 



-69.369 
-44.551 
-2.7575 
155.996 
-705.169 
-696.560 
-0.95653 
-51.7077 



-33.175 
-33.175 
-2.0606 
56.5082 
-611.848 
-611.848 
-0.95456 
-31.5405 



-35.246 
-33.179 
-2.0668 
54.7216 
-612.332 
-611.377 
-0.95538 
-32.7350 



-41.591 
-33.234 
-2.0893 
49.7064 
-613.788 
-609.957 
-0.95785 
-36.3085 



-52.721 
-33.522 
-2.1332 
41.5131 
-616.101 
-607.442 
-0.96219 
-42.3087 



volume energy to about 50 MeV for surface symmetry energy coefficient. The neutron skin size dependence of surface 
symmetry energy is much larger than the skin dependence of volume symmetry energy with A dependence included 
even for a large A of over 200. We can also see that the neutron skin size dependent coefficients for volume energy 
and surface energy Ev(T) s k and Es(T) s k have opposite sign compare to the corresponding skin size independent 
coefficients Ey(T) and E$(T), while volume symmetry energy and surface symmetry energy SV(T) and Ss(T) have 
the same sign for neutron skin size dependent coefficients and independent coefficients. 

The kinetic energy contribution to the neutron skin size dependent and independent coefficients have a much 
more sensitive dependence on the Skyrme parameter set we used. The magnitude of kinetic energy contribution 
follows somewhat the magnitude of the effective mass of the Skyrme parameter set used. Since the neutron skin size 
dependent coefficients Ei(T) s k themselves are somewhat insensitive to the Skyrme parameter set used, the potential 
energy contribution to the coefficients, which are function of density, also are sensitive to the parameter set used. 
The kinetic energy contribution to the neutron skin size dependence Ei(T) s k are now opposite in sign to the kinetic 
energy part to the skin size independent coefficients Ei(T) of energy expansion for all the terms. In turn, they 
are all opposite sign to the skin size independent coefficients Ei{T) themselves. The kinetic energy contribution to 
the neutron skin size dependent volume symmetry energy coefficient Sv{T) s k has a small magnitude similar to the 
magnitude of the kinetic energy contribution to the skin size dependent volume energy coefficient Ev{T) s k- The 



kinetic energy contribution to the surface symmetry energy coefficients, both the neutron skin size dependent and 
independent ones Ss(T) s k and Sg{T), have largest magnitude among the energy expansion coefficients. It is much 
larger even with A dependence included than skin size dependence of volume symmetry energy Sv(T) s k and other 
energy expansion coefficients. 

When the Helmholtz free energy is minimized, from the values of "T-indp" and "T 2 " for T = in Table U the 
temperature and neutron skin size dependence of the energy at low T becomes 



E(A,Z,T,t) 



for SLy4 parameter set, 

E(A,Z,T,t) 



for SkM* parameter set, and 
E(A,Z,T,t) 



(15.308 - 0.012T 2 )A + (20.008 + 0.548T 2 )A 2/3 

+ (31.113 + 0M5T 2 )I 2 A - (41.035 + 2.577T 2 )/ 2 ^ 2/3 
(5.520 - 0.007T 2 )yl - (17.077 - 0.020T 2 )A 2/3 
+ (63.203 + 0.031T 2 )/ 2 A - (711.559 + 0.916T 2 )/ 2 v4 2/3 

+E C (T, t)-^ + E dlf (T, t) — + E ex (T, t) 



(15.127 - 0.019T 2 )A + (18.756 + 0.537T 2 )A 2/3 

+ (29.655 + 0.548T 2 )/ 2 A - (43.596 + 2.560T 2 )/ 2 A 2/3 
(5.584 - 0.008T 2 )A - (17.261 - 0.021T 2 )A 2/3 

+ (62.787 + 0.030T 2 )/ 2 A - (706.246 + O.MOT 2 )I 2 A 2/3 
Z 2 Z 2 z 4 /' 3 

+E C (T, 4)^3+ E dlf (T, t)— + E ex (T, t) 



-(15.310 - 0.039T 2 )A + (18.303 + 0.500T 2 )yl 2/3 

+ (19.685 + 0.423T 2 )/ 2 A - (33.175 + 2.061T 2 )/ 2 A 2/3 
(6.981 - 0.008T 2 ),4 - (21.691 - 0.022T 2 )A 2/3 

+ (63.366 + 0.031T 2 )/ 2 A - (611.848 + 0.955T 2 )/ 2 ^ 2/3 
Z 2 Z 2 z 4 i 3 

+E C (T, t)j m + E diJ (T, t)— + E ex (T, t) 



A 1 / 3 



(8) 



A 1 / 3 



(9) 



A 1 / 3 



(10) 



for SkM(m* = m) parameter set in MeV and fm units. 

From the values of "T-indp" and "T 2 " for T = in Table HI the temperature and neutron skin size dependence of 
the kinetic energy at low T becomes 

E K (A,Z,T,t) = [(28.653 + 0.012T 2 )A - (32.106 - 0.548T 2 )A 2/3 

-(24.674 - 0.545T 2 )/ 2 A + (171.020 - 2.577T 2 )/ 2 A 2/3 

-(3.878 + 0.007T 2 ),4 + (12.325 + 0.020T 2 )A 2/3 

+ (1.332 + 0.031T 2 )/ 2 A - (276.281 + 0.916T 2 )/ 2 A 2/3 

for SLy4 parameter set, 

E K (A,Z,T,t) = [(25.483 + 0.019T 2 )A - (26.275 - 0.537T 2 )A 2/3 



A 1 / 3 



(11) 
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-(33.636 - 0.548T 2 )/ 2 ^ + (166.818 - 2.560T 2 )/ 2 A 2/3 
-(0.489 + 0.008T 2 )A + (1.409 + 0.021T 2 )v4 2 / 3 
-(1.340 - 0.030T 2 )/ 2 A - (40.353 + 0.940T 2 )/ 2 A 2/3 



A 1 / 3 



(12) 



for SkM* parameter set, and 

E K (A, Z,T,t) 



(19.957 + 0.039T 2 ) A - (15.845 - 0.500T 2 )A 2/3 

-(6.881 - 0.423T 2 )/ 2 yl + (56.508 - 2.061T 2 )/ 2 A 2/3 
(0.387 + 0.008T 2 )A + (1.123 + 0.022T 2 )A 2/3 
-(0.871 - 0.031T 2 )/ 2 yl - (31.541 + 0.955T 2 )/ 2 A 2/3 



A 1 / 3 



(13) 



for SkM(m* = m) parameter set. 

Tableland above Eqs.([51)- (fT3")) show that the contribution of the T 2 dependent term in kinetic energy (see Eq.(|X3|) 
to the neutron skin size dependent coefficients Ei(T) s k are very small which is consistent with the result in Ref.[l3j. 
In Ref . fl3| , we estimated approximately the neutron skin size dependence coming from the T 2 term in kinetic energy 
(see Eq. (|A3p ) and found the effects were small. The contribution of the T 2 dependent term in kinetic energy to 
the neutron skin size are insensitive to the Skyrme parameter set used and to the temperature. The neutron skin 
dependent volume energy coefficient Ev(T) s k has smallest effect of T 2 dependence and the surface symmetry energy 
Ss(T) s k has a largest effect similar to the T 2 dependence of the corresponding skin independent coefficients of energy 
expansion. 

By fitting the values of Ei(T) s k for temperatures which are T = 0, 1, 2, and 3 MeV in Table HI the T-dependences 
of the neutron skin dependent coefficients are 



5.52016 - 0.00722T 2 ] 

5.52016 - 0.00008T - 0.01005T 2 - 0.00007T 3 
-17.07720 + 0.01959T 2 ! 



Ev(T) s k 
E s (T) sk 

S v {T) sk [= 63.20304 + 0.03095T : 

S s {T) sk [= -711.5588 - 0.91606T 2 ] 



= -17.07720 + 0.00091T + 0.035997^ + 0.00044T 3 



= 63.20304 + 0.00560T + 0.01927T 2 - 0.00322T 3 



-711.5588 + 0.00963T - 0.01940T 2 + 0.06117T 3 



(14) 
(15) 
(16) 
(17) 



for SLy4 parameter set. The first expressions given in the square parenthesis are from the values for T = in Table 
Usame as in Eqs.(|5))- ([1"3")) for comparison. 



E v {T) sk [= 5.58368 - 0.00761T 2 ] 

E s (T) sk 
Sv(T) s k 

S s {T) sk [= -706.2461 - 0.94030T 2 ] 



5.58368 - 0.00042T - 0.01018T 2 - 0.00013T 3 
-17.26147+ 0.02082T 2 ] 

-17.26147 + 0.00154T + 0.03763T 2 + 0.00056T 3 
62.78658 + 0.02954T 2 ] 

62.78658 + 0.00634T + 0.04623T 2 - 0.00358T 3 



= -706.2461 + 0.06100T - 0.15020T 2 + 0.08320T 3 



for SkM* parameter set. 



E v {T)sk [= 6.98087 - 0.00793T 2 
E s {T) sk [= -21.69093 + 0.022014T 2 



= 6.98087 - 0.000120T - 0.01048T 2 - 0.00008T 3 



-21.69093 + 0.00070T + 0.03447T 2 + 0.00026T 3 



(18) 
(19) 
(20) 
(21) 

(22) 
(23) 



S v (T) sk [= 63.36637 + 0.03144T 2 ] 

= 63.36637 - 0.00421T + 0.02050T 2 - 0.00150667T 3 
S s (T) sk [= -611.8476 - 0.95456T 2 ] 

= -611.8476 + 0.03905T - 0.54265T 2 + 0.01900T 3 

for SkM(m* — m) parameter set. 

By fitting the values for T = 0, 1, 2, and 3 MeV in Table HI the temperature T and neutron skin 
of the expansion coefficients of the kinetic energy are 



E v (T,t) 



K 



E s (T,t) K 



Sv(T,t) K 



Ss(T,t) K 



for SLy4 parameter set. 

E v (T,t) K 



Ss(T,t) 



K 



for SkM* parameter set. 

E v {T,t) K 



-32.10626 + 0.01492T + 0.45071T 2 + 0.00230T 3 
+ (12.32541 - 0.00188T - 0.03149T 2 - 0.00207T 3 

-24.67405 + 0.54454T 2 ) + (1.33205 + 0.03095T 2 
-24.67405 - 0.13177T + 0.26686T 2 - 0.02916T 3 ) 

-(1.33205 + 0.03334T - 0.00683T 2 - 0.00072T 3 )-4-- 

A 1 ' 3 

171.0202 - 2.57693T 2 ) + (-276.2811 - 0.91606T 2 ; 



171.0202 + 1.18275T - 2.36020T 2 + 0.23595T 3 ) 
+ (-276.2811 - 0.18403T + 1.15055T 2 + 0.04358T 3 



Es(T,t) K [= 



S v (T,t) K [= 



AV3 J 



28.65250 + 0.01151T 2 ) + (-3.87804 - 0.00722T 2n ' 
28.65250 - 0.00801T - 0.04819T 2 - 0.00264T 3 ) 
F(-3.87804 + 0.00117T+ 0.01064T 2 + 0.00078T 3 )— ^ 

-32.10626 + 0.54833T 2 ) + (12.32541 + 0.01959T 2 ^ * 



A 1 / 3 



A 1 / 3 ' 



t 



A 1 / 3 ' 
t 

] A~U 3 ~ 



25.48323 + 0.01925T 2 ) + (-0.48934 - 0.00761T 2 )— 

A 1 ' 3 



25.48323 + 0.00253T - 0.04169T Z - 0.00038T J 



-(-0.48934 + 0.00007T - 0.01019T 2 - 0.00003T 3 )— 

-26.27463 + 0.53708T 2 ) + (1.40926 + 0.02082T 2 )— r^l 
' K A 1 ' 3 

-26.27453 - 0.00990T + 0.47594T 2 - 0.00145T 3 ) 

-(1.40926 - 0.00022T+ 0.02858T 2 + 0.00008T 3 )— ^ 

v A 1 / 3 

-33.63624 + 0.54816T 2 ) + (-1.33997 + 0.02954T 2 )-^] 
-33.63624 - 0.08819T + 0.24295T 2 + 0.00042T 3 ) 
+ (-1.33997 - 0.00176T + 0.05211T 2 + 0.00159T 3 ) — " T 

166.8184 - 2.55977T 2 ) + (-40.35256 - 0.94030T 2 ) 
166.8184 + 0.07982T - 0.93650T 2 - 0.09752T 3 ) 
-(-40.35256 - 0.00212T - 1.22403T 2 - 0.01232T 3 



t 

A^ 3 
t 



A 1 / 3 ' 
t 



A 1 / 3 



(19.95733 + 0.03948T 2 ) + (-0.38673 - 0.00793r 2 )-^-^] 
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(19.95733 + 0.00161T + 0.01917T 2 - 0.00003T 3 ) 



t 



E s (T,t) K 



+ (-0.38673 + 0.00006T - 0.00962T 2 - 0.00002T 3 )— - 

A 1 ' 3 

(-15.84463 + 0.49952T 2 ) + (1.12277 + 0.02201T 2 ) ' 1 



(34) 



-15.84463 - 0.00388T + 0.46826T 2 + 0.00162T 3 
(1.12277 - 0.00024T + 0.02672T 2 + 0.00005T 3 



A 1 / 3 ' 
t 



A 1 / 3 



(35) 



Sv(T,t) 



K 



Ss(T,t) K 



(-6.88128 + 0.42271T 2 ) + (-0.87085 + 0.03144T 2 )— ^] 
(-6.88128 + 0.01595T + 0.34709T 2 - 0.00140T 3 ) 

+ (-0.87085 + 0.00273T + 0.04486T 2 + 0.00089T 3 )— ^- 

A 1 ' 3 

t 



(36) 



(56.50820 - 2.06058T 2 ) + (-31.54054 - 0.95456T 2 ) 
(56.50820 - 0.15544T - 1.63960T 2 + 0.00843T 3 ) 
+ (-31.54054 - 0.02086T - 1.16564T 2 - 0.00796T 3 



A 1 / 3 
t 

) aU 3 ~ 



(37) 



for SkM(m* — m) parameter set. 

Comparing the first line (or Eqs.(|8"|)- flT3"|) ') and second line of above Eqs. (fT4")) - (|3"7|) . we can see the temperature T 
dependence of the neutron skin size dependent coefficients Ei(T) s k comes not only from the T 2 term in the kinetic 
energy, Eq. (|A3j) . but also from the potential energy through the different saturation density for different temperature. 
Especially, the neutron skin size dependence of kinetic energy expansion coefficients for SLy4 parameter set have 
opposite sign in their T dependence compared with the T 2 term in kinetic energy (T 2 dependence in Eqs.(15|- (|13[) ). 
Table|I]and Eqs. (Ti"4)) -([2"5 j) show that the T dependence of the neutron skin size dependent coefficient is much faster than 
the T 2 term in kinetic energy of Eq. (|A3[) alone for volume energy coefficient Ey(T) s ). and surface energy coefficient 
Es(T) s k- By contrast this it is much slower than the T 2 term for volume symmetry energy coefficient Sv(T) s k and 
surface symmetry energy coefficient Ss{T) s k except for volume symmetry energy Sv(T) s k of the SkM* parameter set. 

Eqs. (j2"6"l) - ([3"7"1) show that even the kinetic energy also has an extra T dependence, beside the T 2 term in kinetic energy 
of Eq. (j A3 jl , through the T dependence of saturation density. For SLy4 Skyrme interaction, the T dependence of the 
kinetic energy part of the neutron skin size dependent energy expansion coefficients Ei(T) s k has opposite sign with 
the T 2 term of kinetic energy Eq. (|A3|) . For this interaction the T dependence of the kinetic energy part of the volume 
symmetry energy coefficient Sy(T) s k is rather linear as compared to a T 2 behavior. For SkM* and SkM(m* = m) 
interactions, the kinetic energy part of the neutron skin dependent energy coefficients E i (T) s i : has the same sign in its 
T dependences with the T 2 term in the kinetic energy of Eq. (IA3|) but has a faster dependence of T compared to the 
T 2 term of Eq. (|A3|) . The neutron skin independent kinetic energy expansion coefficients have a slower T dependence 
than a T 2 term for a kinetic energy, Eq. (|A3p . except the volume energy coefficient Ey(T, t — 0)k for SLy4 and SkM* 
parameter sets which have an opposite sign compared to T 2 term of kinetic energy Eq. (|A3[) . The neutron skin size 
independent volume symmetry energy coefficient Sy(T,t = 0)k for SLy4 parameter set has a linear T dependence 
comparable order to the T 2 dependence. 

The kinetic energy expansion at zero T, from Eqs. pT|) - P^)l or from Eqs. (|2"6"|) - (f3~T|) , are 



E K (A,Z,T = 0,t) = 

for SLy4 parameter set, 

E K (A,Z,T = 0,t) = 

for SkM* parameter set, and 

E K (A,Z,T = 0,t) = 



28.653A - 32.106^ 2/3 - 24.674/ 2 ^l + 171.020/ 2 A 2/3 
3.878A+ 12.325A 2/3 + l.332I 2 A - 276.281/ 2 yl 2/3 

25.483,4 - 26.275A 2/3 - 33.636/ 2 ^ + 166.818/ 2 A 2/3 
+ [-0.489A + 1.409A 2/3 - 1.340/ 2 A - 40.353/ 2 A 2/3 



19.957A - 15.845A 2/3 - 6.881/ 2 A + 56.508/ 2 A 2/3 
0.387 A + 1.123A 2 / 3 - 0.871/ 2 A - 31.541/ 2 A 2 / 3 



A 1 / 3 



(38) 



A 1 / 3 



(39) 



A 1 / 3 



(40) 
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for SkM(m* = m) parameter set. For an infinite nuclear matter, the surface terms disappear and only the volume 
terms survive. These results show that the surface symmetry energy coefficient has the largest effect from kinetic 
energy as compared to the other coefficients. Here we can see the surface kinetic energy and the volume symmetry 
kinetic energy coefficients are negative. However the total kinetic energy and total symmetry kinetic energy including 
A factors are positive. The neutron skin dependent kinetic energies with A factor included are negative. Compre to 
result from Fermi gas model, 

E K = 12I 2 A + 9I 2 A 2/3 (41) 

which is good for high T or low density limit without any interaction. Here both the volume and surface symmetry 
energies are positive. With Skyrme interaction, the isospin dependent part of the effective mass in a finite nuclei 
may become negative depending on the force parameter and densities of proton and neutron. Thus the signs in 
Eqs.flHHJ-lHU result - 

IV. CONCLUSION AND SUMMARY 

Understanding properties of the symmetry energy is important in many area of nuclear physics as mentioned in 
the introduction. In this paper we studied properties of the symmetry energy, both volume and surface parts, along 
with other terms which appear in the Weizsacker mass formulae. Our investigation was based on a finite temperature 
density functional approach. In a finite temperature approach, the dependence of various terms on temperature can 
be obtained. Energetic probes lead to excited nuclei which may be characterized by a hot liquid drop extension of the 
Weizsacker mass formulae. We used several different interactions of the Skyrme type to examine the dependence of 
various quantities on the interaction and associated effective masses that appear. Our analysis in the present study 
emphasized the role of the neutron skin on various terms that appear in the mass formulae. The radii of protons and 
neutrons were therefore allowed to be different in a Saxon- Wood form for the density distributions of these particles. 
We then proceeded to calculate various terms using an expansion about the equal radii point. The corrections that 
arise from a neutron skin are then proportional to the skin thickness t over the radius R or t/R ~ t/A 1 ^ 3 . The skin 
thickness t/R itself can be proportional the neutron excess fraction I = (N — Z)/A when the proton and neutron 
central densities are the same. Thus the A and / dependences of various terms in the mass formulae are modified. 

Table U contains the results for three Skyrme interactions, SLy4, SkM* and SkM(m* = m). Results for the volume 
energy Ey{T) and Ev{T) s k, surface energy Es(T) and Es(T) s k, volume symmetry energy Sy(T) and Sv(T) s k, and 
surface symmetry energy Ss{T) and Ss(T) s k are given. The terms with an additional subscript "sfc" are the skin 
coefficients of Eq.©. The kinetic energy contributions, labeled "Kine", to each term are also given. The difference 
of the total and kinetic term is from the interaction. 

The results show that the neutron skin size dependent and independent energy expansion coefficients are rather 
insensitive to the Skryme interaction used while the kinetic energy and potential energy expansion coefficients sep- 
arately are sensitive to the interaction used and somewhat follow the effective mass of the Skyrme parameter. The 
temperature dependence of the neutron skin size dependent energy expansion coefficients are much more insensitive 
than the temperature dependence of neutron skin size independent coefficients. The magnitude of the neutron skin 
size dependent coefficients are largely different for different coefficients compared to the neutron skin size independent 
coefficients. The neutron skin size dependent volume energy coefficient has the smallest magnitute while the surface 
symmetry energy coefficient has the largest magnitude. The neutron skin size dependence of surface symmetry energy 
is much larger than the skin size dependence of the volume symmetry energy with A dependence included. The surface 
symmetry kinetic energy coefficients, both the neutron skin size dependent and independent ones, have the largest 
magnitude of the kinetic energy expansion coefficients. The neutron skin size dependent volume energy coefficient has 
the smallest temperatutre dependence and the neutron skin size dependent surface symmetry energy coefficient has 
the largest temperature dependenc similar to the neutron skin independent coefficients. The temperature dependence 
of the neutron skin size dependent coefficients is smaller than the temparature dependence of the skin size independent 
coefficients. The temperature dependences of the neutron skin size dependent energy coefficients have a large effect 
from the different saturation density for different temperature and thus do not follow T 2 behavior of the explicit T 2 
dependence of kinetic energy. 

Considering the neutron skin size dependence t/R factor with R = r^A 1 / 3 , the neutron skin dependent energy 
expansion coefficients has an extra A -1 / 3 factor in the energy expansion. On the other hand, if we relate the 
dimensionless factor t/R to the isospin factor / = (N — Z)/A then the energy expansion including the neutron skin 
size dependence introduce an extra / factor. With this extra / dependence we may be able to extract some information 
on the neutron skin size by expanding the empirical energy of various nuclei with including odd power of / up to third 
order. 
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Appendix A: Skyrme interaction 

The Hamiltonian for a Skyrme interaction is 



H(r) = H B {r) + H s (r)+H c {r) 



Hu = 



H s = 



2m 3 
1 

+ 4 



*o 
2 

h 

12 

1 r 

16 

1 

~16 

D 2 



-Tp + 



h z 



+ 



2m„ 

(l + y) P 2 -(l + *o) (pI + pI) 

X i + x i)p 2 -{\+A {pI+pD 



1 








pT ~4 


*i 




-* 2 Q+* 2 ) 



(PpTp + p n T n ) 



H c = -ttPp 



3tl ( 1 + |)_ (2 ( 1 + |)] ( vp, 

[3tl(1+f) _ (2(1+f) 

(r) J d 3 r 



1 

16 



p\7 2 p + 



16 



3*i 
3*i 



Xi + * 2 



zi + *2 



X2 



X2 



[{Vppf + (V P „) 2 ] 

(p p V 2 p p + Pn^Pn) 



,3 J Pv^f 



3e^ / 3 
4 Itt 



P 4 p /3 (r) 



(Al) 



2 ' ' ' ' j \r- f' 

The H{r) has a bulk part H B (r), a surface part Hs{r) with gradient terms and a Coulomb term H c {r). The effective 



mass to" is 



m 2m [ 1 

— = 1 + T2 



h (l + s) +fa (l + 2 



1 

p -l 



h I ^ +a;i j -* 2 Q + .x 2 



P 9 



2m f 1 i 
1 + ~W 1 16 [3tl + (5 + AX2)h] PTz 



p(2y - 1) 



where q = n,p for neutron or proton. At low T, 



3 /6tt' 

7 



2\ V3 



pT 



5ir 2 m* 2 / 7 n4/3 
3^ 4 Vfo^J 



P y 3 T 2 + . . . 



(A2) 



(A3) 



The gradient terms are important in finite nuclei and the Coulomb term is important for the charged proton compo- 
nent. The * 0j *i) *2, *3 and x , x\, Xi, x 3 are parameters. Different choices of these parameters give rise to different 
Skyrme interactions. 



Appendix B: Expansion of density 



The central density p qc of the density distribution for a given value of R q should be determined to give a fixed 
number of nucleons N q ; 



N q (R q ) = J d 3 rp q (r)=4ir J r 2 dr- 



PqcjRg) _ fJf A R 3 
+ e (r-R q )/a -PvcyKq) 3 H q 



l + TT 2 



This normaization condition gives up to first order in dR q = (R q — R), 



dR q + p qc -^-R 2 







3 + 7T 2 









dR q = 



(Bl) 



(B2) 
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dp qc 


Pqc(R) 


3+**(a) a ; 


dR 


R 


"i+- 2 (t) 2 " 



Pqc(Rq) = Pqc{R) 



1 - 



3W(f) 2 



(B3) 
(B4) 
This 



Thus the first order correction coefficient to p™ due to fixed JV g is the zeroth order term times — ^ 

result is independent of which type of particle, neutron or proton. 

If we keep the particle number N q and the total central density p c — p nc + p pc to be a constant while varying R q , 
then 



p c (R) = Pnc(Rn) + Ppc(Rp) 
= p n c(R) + Ppc(R) - 
= Pnc(R) + Ppc(R) 



3 + ^ 2 (f)^ 



Pnc(R) 



Rn — R 

R 



Thus we get 



and 



R n — R 
R 



Ppc 



Rp — R 
R 



= 



Rp — R 



t = R n 



Rv 



——(Rn - R) 
Ppc 

{R n — R) — {R p -R) = (R n -R)(l 



Pnc 
Ppc 



Finally we have 



R 



Ppc, 



V 



R = 



Pnc 



Ppc(R) 



Rp — R 
R 



Ppc 



-{Rn — R) 



(B5) 



(B6) 



(B7) 



(B8) 



to lowest order in t. The same result can be obtained by requiring A = iV + Z constant with keeping the central 
densities p qc unchanged in the dR q expansion. When R q = R, p qc (R) / p c (R) — N q /A and thus we get 



t-n Rn R j^-, 



tp — Rp 



N 

R = — -t 
A 



(B9) 



This is the same result given in Ref.|2l]. 

Since the size R q depends on the central density p qc for a given value of particle number N q as in Eq. (|Bip . the 
neutron skin size t is related to the particle number N q and the central density p qc - Using Eqs. (|Bl|) and (|B8|) we can 
obtain following conditions to lowest order in t. 

Air 

N-Z = —[(R 3 n + 7T 2 a 2 Rn)pnc~(Rl + TT 2 a 2 Rp)pp C ] 
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t 

R 



-ir z a z R 1 + 



Ppc ^ 



in 



ir 2 a z R 1 



in 



Pc R 



Ppc 



3 

47T , 

T 



.4 

1 1 

3 'A 



Pnc Ppc 



3 

Pc 



Pnc 



^- (R :i + vrVi?) ( Pnc - Ppc ) + ^- (3i? a + n z a z R) ( ^ Ppc + ^p nc ) - 



R 



in 

T 



o i Pnc Ppc 
«J I Ppc T" Pnc 

V Pc Pc 



A ( Pnc Ppc 
*-* I Ppc ~T~ Pnc 

Pc \ Pc Pc 



2 1 



7ra i 
v R I 



t 

R 



Pc 



1 

3 V 2 (! ~ 2/c)2/c 



Pnc 

N - Z 



A 



(N-Z)-A 
(1 - 2y, 



Ppc 



l+(f) 2 
1+1(1)' 



1 + (1) 



3 U. 



where y c = p pc /p c is the proton fraction of the central density. For the case of p nc ~ Ppc, Uc 
(1 — y c ) = 1/2 — e. Then the factor [2(1 — y c )y c ]~ 1 becomes 



t 

R 



(BIO) 



(Bll) 
1/2 + e and 



2(1 - Vc)Vc 



(l-2e)(l + 2e) 



2(1 -4e 2 



(B12) 



Thus for an uniform distribution (diffuseness parameter a = 0) with p nc w /9 pc , the neutron skin size t/R of Eg. dBlip 
becomes, up to first order in e, 



t 

R 



2N 



3 A -3 (1 ^ c) 
which is the result given in Ref.[21]. On the other hand, as another form, 



(B13) 



N 
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— — = 4f [(Rl + 7rVi?„) - (R 3 P + n 2 a 2 R p )] 

Pp 
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.4 
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Pc P 
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~p7# 



TrVi? | 1 
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7c~R 
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7c~R 



-ir 2 a 2 R\ 1 



ir 2 a 2 R | 1 



Pnc t 

Pc R 

Pnc £ 



n 2 a 2 R 



f Ppc _ 




t 


4tt 


R 3 H 


" 3 R 


t 


V Pc 


Pc /_ 


i? 


" Y 6 






R 



t 

R 



1 



N 



rra \ 
, R I 



t 

R 



1 Pc 
3 A \p n 



z 

Ppc 



Tra) 
, fl J 



3 \ R 



1 1 

34 V 1 



N 



Z 

Vc 



, R 



,2 1 



1 + 1(1) J 



(B14) 
(B15) 



Even if Eqs. (|Bll[) and (|B15|) look different they are the same equation. The proton ratio y c has the range of 
Z/A < y c < 1/2 for finite nuclei where y c = Z/A when R n = R p and y c — 1/2 for p nc = Ppc- 
For one extreme case of the same central density p nc — Ppc-, the proton ratio y =1/2 and 



t 

R 



2 (N- Z 

3 V A 



i+Hir 



' ira ' 



i + Hf) 2 



(B16) 



Thus the neutron skin size is linearly proportional to the isospin factor / when the proton and neutron central densities 
are same. For the other extreme case of the same size R n = R p = R, the proton ratio is y c = Z/A with 1 — y c = N/A, 
and thus the neutron skin size becomes t/R = 0. However when R n ~ R p with y c = Z/A + e, we have 1 — y c = N/A — e 
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and, from Eq. (|Blip . 



t 

R 



1 



2{N/A-e){Z/A 
1 



TV - Z 
A 



N 
1 



Z 
A 



2{N/A){Z/A) 



A 

1+ N' 



N - Z 
A 



N - Z 
A 



1 + 



7m 
. R 



,2 "1 



1 / A 2 
3 \~NZ 



i+iir 



up to the first order in e. 

Using Eq. (|B4|) . the Fermi density Eq.([T]) is expanded up to first order in (R q 



R) as 
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(B17) 



(B18) 
(B19) 



B20) 



From this result it is easy to show that the quantity p™{r) + p^ir) has a first order correction from skin size t and 
the first order correction vanishes for m = 1. That is the total density p(r) = p n {r) + p P (r) is independent of the 
neutron skin size t = R n — R p up to first order. Thus only the explicit p q dependent terms in Skyrme interaction, not 
the total density p dependent terms, depend on the skin size t up to the first order. 



Appendix C: Integral of density functional 
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Since F[p q ] = J d 3 rf(p q (f f )) for a Feremi density p q (r) 



p q only, can be integrated exactly as a function F(R q ) of R q [13j, we can expand F(R q ) easily in terms of t q — R q — R 
or x q = —t q /a around F(R). That is 



(CI) 
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This is much simpler than using preveous method of expanding in t (Eq. (|B18[l ) first then integrate the results. For 
fixed N q , the central density is expanded as 
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In Ref.[l3j], the Coulomb exchange term is shown only upto 0th order in a/R but the actual calculation included all 
terms of a/R up to 3. For the term with mixed densities, we need to integrate after expansion. Up to 1st order in 
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x = —t/a, 
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For SLy4 parameter with a = 1/6 and ra = 2 
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Here 36/13 = 3 x 2/(13/6). For SkM(m* = m) parameter with a = 1 and m = 2, 
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Here 2 = 3x2/3. For T-independent term in kinetic energy with a = 1 and m — 5/3, 
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Here (15/8) = 3 x (5/3)/(8/3). For kinetic energy we cannot use this method since it uses numerical integration. 

The kinetic energy has term with the form of p™ra* 2 where the effective mass has the form oim q /m* = l + ap + bp q . 
Since 
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and p is t independent up to the first order in the neutron skin t, 
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Since the m* factor depends on p q , the first order term in x = —t/a of (p™m*™ + p™m* n ) does not become zero for 
any m for non-zero n. Similarly total kinetic energy r = t„ + r p has nonzero first order term in the neutron skin t in 
contrast to the total density p = p n + p p which is independent to t up to first order. Total kinetic energy is 
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In Weizacker mass formular it might be better expanding in terms of (R n — R p )/R rather than (R n — R p )/a since 
the first one is independent of A while the second one is dependent on A. 
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